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Abstract 



A search for the light pseudoscalar sgoldstino production in three-body K~ decay K~ — > 7r~7r° P has been 
performed with the "ISTRA+" detector exposed to the 25 GeV negative secondary beam of the U-70 proton 
synchrotron. No signal is seen. Upper limits for the branching ratio Br(K~ — ► ir~TT°P) at 90% confidence level 
vary between 2.0 • 10~ 5 and 0.5 • 10 -5 in the effective mass nip range from up to 190 MeV. Our results 
improve the limits published by the E787 Collaboration in the mass interval between and 120 MeV and are 
the first ones at higher masses. 



1 Introduction 

In supersymmetric models with spontaneous supersymmetry breaking superpartners of a Goldstone fermion 
- goldstino, pseudoscalar P and scalar S, should exist. In various versions of gravity mediated and gauge 
mediated theories one or both of these weakly interacting bosons — sgoldstinos — are light and therefore can 
be observed in kaon decays. Moreover, if sgoldstino interactions with quarks conserve parity,(as in left-right 
extensions of MSSM), and P is lighter than S, so that ms > (tuk — rn n ) and mp < (mx — 2m,,-), sgoldstinos 
show up in the decay K — > ttttP, rather than in much better constrained K — > nS. The phenomenology of light 
sgoldstinos in this scenario is considered in details in [1]. Under assumption that sgoldstino interactions with 
quarks and gluons violate quark flavor and conserve parity, low energy interactions of pseudoscalar sgoldstino 
P with quarks are described by the Lagrangian: 

L = -P ■ (h? ■ dirfdj + hij ■ Titrfuj) (1) 

with 

di = {d,s,b) , m = {u,c,t) , (2) 
with coupling constants hij proportional to the left-right soft terms in the matrix of squared masses of squarks: 

- (LR)2 - (LR)2 

,D _ m D,ij h U _ m U,ij /o\ 



Here the scale of supersymmetry breaking is denoted as v F. The constraints on the flavor-violating coupling 
of sgoldstino to quarks is evaluated using the K° L — K° s mass difference and CP violating parameter e in the 
neutral kaon system: \h® 2 \ < 7- 10 -8 ; \Rehf 2 ■ Imhf 2 \ < 1-5- 10~ 17 . It has been shown in [1] that, depending on 
the phase of sgoldstino-quark coupling, these constraints result in the following limits on the branching ratio: 
Br(K~ — » 7T~7r°P) < 1.5 x 10~ 6 t4x 10 -4 (see Fig. 1 for the diagram of the K~ decay). A search for P in 
charged kaon decays is of particular interest for the case Rehf 2 ~ 0, when the corresponding branching ratio of 
Kl is small. 

Light sgoldstino decays into two photons or into a pair of charged leptons, two photon decay dominating almost 
everywhere in the parameter space. Depending on the parameter g 7 — 2 -^^p a -, where M 77 is the photino 
mass, sgoldstino decays either inside or outside the detector. In the present search we assume that sgoldstino 
decays outside the detector, i.e is invisible. The existing limits on the branching Br(K~ — ► 7r~7r°P) are at the 
level of 4 • 10 -5 [2], whereas the limits for the scalar sgoldstino S can be estimated from the studies of the 
K + — ► -k + vv at the level of 4.7 • 10~ 9 , see [3] for recent review. 

The aim of our present study is to search for invisible pseudoscalar sgoldstino in the decay of K~ . Experi- 
mental setup and event selection are described in Section 2, the results of the analysis are presented in Section 3 
and the conclusions are given in the last Section. 



2 Experimental setup and event selection 
2.1 Experimental setup 

The experiment is performed at the IHEP 70 GeV proton synchrotron U-70. The spectrometer ISTRA+ has 
been described in some details in our recent papers on K e3 [4], [5] and tt~tt°tt° decays [6], here we briefly 
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recall characteristics relevant to our analysis. The 40 m spectrometer ISTRA+ is located in the negative 
unseparated secondary beam line 4A of U-70. The beam momentum is <~ 25 GeV with Ap/p <~ 2%. The 
admixture of K~ in the beam is ~ 3%, the beam intensity is <~ 3 • 10 6 per 1.9 sec U-70 spill. A schematic view 
of the ISTRA+ setup is shown in Fig. 2. The momentum of a beam particle deflected by the magnet Mi is 
measured by four proportional chambers BPCi — BPC4 with 1 mm wire step, the kaon identification is done 
by three threshold Cerenkov counters Ci — C3. The 9 meter long vacuumated decay volume is surrounded by 
eight lead glass rings used to veto low energy photons. The same role is played by 72-cell lead glass calorimeter 
SP2. The decay products deflected in the magnet M 2 with 1 Tm field integral are measured with 2 mm step 
proportional chambers PCi — PC3, with 1 cm cell drift chambers DCi — DC3 and finally with 2 cm diameter drift 
tubes DTi — DT4. A wide aperture threshold Cerenkov counter C4 is filled with He and used to trigger electrons. 
SPi is a 576-cell lead glass calorimeter, followed by HC — a scintillator-iron sampling hadron calorimeter. MH 
is a 11x11 cell scintillating hodoscope, used to solve ambiguity for multitrack events and improve the time 
resolution of the tracking system, MuH is a 7x7 cell muon hodoscope. 

The trigger is provided by scintillation counters Si — S5, beam Cerenkov counters and by the analog sum of 
amplitudes from last dinodes of the SPi : T=Si ■ S2 • S3 ■ S4, ■ Ci • C2 • C3 • S$ ■ E(SPi), here S4 is a scintillation 
counter with a hole to suppress beam halo, S5 is a counter downstream of the setup at the beam focus, X(SPi) 
- a requirement for the analog sum to be larger than a MIP signal. , 

During first run in March- April 2001, 363 million of trigger events were logged on DLT's. During second 
physics run in November-December 2001 350 million trigger events were collected with higher beam intensity 
and somewhat stronger trigger requirements. This information is supported by about 100 million MC events 
generated using Geant3 [7] for dominant K~ decay modes. Signal efficiency for possible sgoldstino production 
has been estimated using one million generated events for the first run configuration and 0.5 million events for 
the second run configuration, for each effective mass mp point in the mass interval from to 200 MeV with a 
step of 10 MeV. These signal events were weighted using the matrix element given in [1]. 

2.2 Event selection 

Data collected in two runs, Spring 2001, first run, and Winter 2001, second run, are used. Some information 
on the data processing and reconstruction procedures is given in [4, 5, 6], here we briefly mention the details, 
relevant for 7r~7r° + missing energy event selection. 

The muon identification (see [5]) is based on the information from the SPi — a 576-cell lead glass calorimeter 
and the HC — a scintillation-iron sampling calorimeter. The electron identification (see [4]) is done using E/p 
ratio — of the energy of the shower associated with charged track and charged track momentum. 

A set of cuts is designed in order to suppress various backgrounds to possible sgoldstino production: 

0) Events with one reconstructed charged track and with two electromagnetic showers in the electromagnetic 
calorimeter SPi are selected. We require the effective mass 111(77) to be within ±40 MeV from m,.. The 
effective two-photon mass distributions for two runs are shown in Fig. 3. Events with vertex inside interval 
400 < z < 1650 cm are selected. 

1) Next cut is a "soft" charged pion identification, tracks having electron or muon flag ( as described in 
[4, 5] ) are rejected. 

2) Events with missing energy Et, eam — E n - — E v o below 3 GeV are rejected. The cut on the missing energy 
serves to reduce K n2 contribution. In Fig. 4 missing energy spectra for the second run data, MC background 
and MC signal with mp of 90 MeV are compared. 

3) Events with 771(77) within ±36 MeV from m,r°135 MeV a selected (Tough ir° cut). 

4) Fourth cut suppresses events with extra photons irradiated by charged particles in a detector material 
upstream M2-magnet. Such photons have nearly the same x coordinate on the SPi face as the charged track, 
i.e event is rejected if at least one photon has \x c h — x 7 | < 10 cm. 

5) Fifth cut removes events having good K e3 2C-fit. 

6) The following cut is a restriction on the charged pion momentum p* in the kaon rest frame: p*(ir~) < 
180 MeV. It is designed to suppress the tails of the K n 2 decay. 

7) Seventh cut is against remaining muon contamination.lt is required that the energy in SPi, associated with 
charged track, should exceed 0.7 GeV. This cut requires "early" hadron shower in SPi and has , unfortunately, 
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low signal efficiency. We plan to get rid of it in our final analysis. 

8) The following cut is designed to suppress 7r~7r°7 contribution: only the events with missing momentum 
crossing SPi transverse plane within SPi working aperture are selected. 

9) The "Veto" cut uses information from the Guard System GS and the guard electromagnetic calorimeter 
SP2: absence of signals above noise threshold is required. The influence of the "Veto" cut on the two-photon 
effective mass spectrum is illustrated in Fig. 5. No ir° signal is seen after the "Veto" cut. 

10) As a result of the previous cuts, 7r°-signal practically disappears from the 77 mass spectrum. This allows 
to perform effective linear background subtraction procedure. The interval ± 18 MeV around m(7r°) mass value 
is used as the signal region, remaining events are used for the background estimate. 

The data reduction information is given in Table 1 for two runs and is compared in Tables 2,3 with MC- 
background statistics and with MC-signal statistics for the sgoldstino mass of m = 90 MeV. 

The influence of the last eight cuts on the missing mass squared spectrum (Pk — P v - — Ptt°) 2 , which is used 
as the signal "estimator" is is shown in Fig. 6 for the second run data. The left wide bump in Fig. 6 is due to 
K^s decays, the shift to negative missing mass squared is caused by the use of the pion mass in its calculations. 
Second peak is caused by 7r~7r°(7r°) decay with gammas from second tt° escaping detection. 

Fig. 7 illustrates the background subtraction procedure: in the first (third) figure the missing mass squared 
spectrum after cut 10 is shown for the first (second) run (solid curve) together with the normalized back- 
ground(dashed curve) estimated from the 771(77) spectrum side bands ( 99-117 MeV and 153-171 MeV). The 
second(forth) figure show signal-background difference for two runs. No sgoldstino signal is observed, the peak 
around 0.02 GeV 2 is due to remaining 7r~7r°(7r°) background. 



Table 1: Event reduction statistics for two runs. 



Cut 


Events run 1 


N l /N l+1 


Events run 2 


N l /N l+1 


(0) 1 7r~ and 111(77) ncar ni(7r°) 


4398982 




11339975 




(1) no (c, fi) 


978752 


1.48 


8762220 


1.29 


(2) E mjs > 3.0 GeV 


716830 


4.16 


1430781 


6.12 


(3) narrow 111(77) band 


581653 


1.23 


1297087 


1.10 


(4) conv. gammas 


322073 


1.81 


723248 


1.79 


(5) no K e3 fit 


230424 


1.40 


487246 


1.48 


(6) p*(n-) < 180 MeV 


150926 


1.53 


322298 


1.51 


(7)E SPl > 700 MeV 


17101 


8.83 


51624 


6.24 


(8) 10 < rr < 60 cm 


12922 


1.32 


40356 


1.28 


(9) Veto 


5255 


2.46 


11183 


3.61 


(10) ±18 mcV 


2691 


1.95 


6125 


1.94 



3 Analysis and results 

In order to calculate the upper limits we have fitted the residual missing mass distribution by the sum of the 
background and the signal Gaussian with fixed width of a = 0.003 GeV 2 , as determined from the signal MC. 

The signal MC is also used to determine the efficiency of our cuts. The background has been described by 
two components: the Gaussian for 7r~7r (7r ) peak and a linear background. For the second run an additional 
wide Gaussian centered at -0.015 has been used. Fit results are illustrated in Fig. 7 for the first and second 
run, respectively. 

The upper limits at the 90% confidence level are calculated from the value 

N UL = N obs + 1.64 • a , (4) 
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Table 2: Event reduction statistics for MC data, first run. 



Cut 


Events mix MC 


Ni/N i+1 


MC signal with m=90 MeV 


N t /N l+1 


(0) 1 7T~ and m(77) near m(7r°) 


1278382 




207315 




(1) no (c, fj) 


1156685 


1.11 


203841 


1.02 


(2) E mls > 3.0 GcV 


158181 


7.31 


168179 


1.21 


(3) narrow 111(77) band 


149679 


1.06 


165105 


1.02 


(4) conv. gammas 


91251 


1.64 


94503 


1.75 


(5) no K e3 fit 


48067 


1.90 


68325 


1.38 


(6) p*(ir-) < 180 MeV 


32151 


1.50 


67528 


1.01 


(7)E SPl > 700 MeV 


6880 


4.67 


22346 


3.02 


(8) 10 < rr < 60 cm 


5283 


1.30 


17883 


1.25 


(9) Veto 


2257 


2.34 


17770 


1.01 


(10) ±18 MeV 


1521 


1.48 


15737 


1.13 



where N a i, s is the number of events in the signal Gaussian. During the fit we impose the constraint that N b s 
is nonncgative. Sigma is a non-Gaussian error in N a b s . The upper limit UL is calculated as 

TT , _ N UL • 0.2116 • 0.98798 • e(K n2 ) 

UL ~ N{K^)-e (5) 

with 0.2116 equal to Br(K n2 ) ; with 0.98798 equal to Br(w° — > 77) and N(K 7r2 ) equal to the number of 
reconstructed ir~ir° decays found to be 1500000 events for the first run and 4500000 events for the second run. 

The values e and s(K n2 ) are respective acceptances for K~ — > 7r~7r° P and K n2 decays. The acceptance 
s(K n2 ) is equal to 17.49 % for the first run and 24.92 % for the second run. This acceptance includes both 
the reconstruction efficiency and geometrical acceptance. The results of the fits are given in Table 4. As the 
background conditions for the two runs are different, the first run provides better upper limits for some mp 
regions regardless of the lower luminosity. The estimated systematic errors for the upper limits are at the level 
of 20 %. 

The final result obtained using combined statistics of two runs is given in Table 5. The weighted ratio 
s(K^ 2 )/e with the weights proportional to the runs statistics has been used for the combined data sample. 
Fig. 9 shows a comparison of our result with that published by the E787 collaboration [2]. 

Special attention should be paid to the region near m^o mass. If we assume that the sgoldstino mass is 
exactly equal to that of tt° and that the remaining tt° peak is due to sgoldstino signal, then the upper limit for 
this point is 2.8 • 10~ 4 . 



4 Summary and conclusions 

A search for a possible pseudoscalar sgoldstino production in the decay K~ — > tt~tt°P has been performed. It 
was assumed that sgoldstino decays outside the setup, i.e is invisible. No signal is seen in the mp mass interval 
between and 200 MeV. The 90% confidence upper limits obtained vary between 2.0 • 10~ 5 and 0.5 • 10~ 5 for 
the sgoldstino mass range of 4- 190 MeV. These results improve the confidence limits published by the E787 
Collaboration. 

In future, we plan to extend our search for a sgoldstino signal in the scenario when the decay P — > 77 
happens inside the setup. 

The authors would like to thank D.S. Gorbunov, V.A. Matveev, V.A. Rubakov, for numerous discussions. 
The INR part of the collaboration is supported by the RFFI fund contract N00-02-16074. 
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Table 3: Event reduction statistics for MC data, second run. 



Cut 


Events mix MC 


Ni/N i+1 


MC signal with m=90 McV 


N t /N l+l 


(0) 1 7T~ and m(77) near m(7r°) 


5117109 




120594 




(1) no (c, n) 


4627156 


1.11 


119292 


1.01 


(2) E mls > 3.0 GcV 


618446 


7.48 


98098 


1.22 


(3) narrow 111(77) band 


565772 


1.09 


96504 


1.02 


(4) conv. gammas 


346226 


1.63 


54206 


1.78 


(5) no K e3 fit 


118536 


2.92 


37124 


1.46 


(6) p*(n~) < 180 McV 


114160 


1.04 


36835 


1.01 


(7)E SPl > 700 McV 


21824 


5.23 


11800 


3.12 


(8) 10 < rr < 60 cm 


16597 


1.28 


9462 


1.25 


(9) Veto 


7229 


2.30 


9407 


1.01 


(10) ±18 MeV 


5114 


1.41 


8541 


1.10 



Table 4: Mass dependence of upper limits UL calculated using two runs separately, left part — first run, right 
part second run. 



M, MeV 


Eff. e, % 


N obs ± <J 


UL 


Eff. e, % 


N obs ± cr 


UL 





0.64 


0±4.80 


1.84 • IO" 5 


0.85 


0±17.51 


2.39 • 10~ 5 


10 


0.66 


0±4.76 


1.76 • 10" b 


0.86 


0±17.90 


2.40 • 10~ b 


20 


0.67 


0±4.65 


1.69 • IO" 5 


0.92 


0±16.81 


2.10 • 10~ 5 


30 


0.74 


0±4.49 


1.48 • 10" b 


0.99 


0±14.88 


1.73 • 10~ b 


40 


0.81 


0±4.32 


1.30 • IO" 5 


1.05 


0±12.32 


1.65 • 10~ b 


50 


0.88 


0±4.19 


1.16 • 10" b 


1.18 


0±9.73 


1.56 • 10~ b 


60 


1.05 


0±4.19 


0.98 • 10" & 


1.32 


0±7.67 


1.10 • IO" 5 


70 


1.09 


0±4.51 


1.01 • IO" 5 


1.73 


0±6.35 


0.70 • 10~ b 


80 


1.21 


0±5.71 


1.15 • IO" 5 


1.62 


0±5.62 


0.66 • 10" 5 


90 


1.33 


0±10.97 


2.02 • IO" 5 


1.76 


0±5.26 


0.57 • 10~ 5 


100 


1.41 


0±24.59 


4.24 • IO" 5 


1.90 


0±5.40 


0.54 • 10~ 5 


110 


1.53 


10.96±18.16 


3.96 • 10" & 


2.09 


0±6.86 


0.62 • IO" 5 


120 


1.60 


17.22±15.99 


2.24 • 10~ 5 


2.27 


0±12.44 


1.04 • 10~ 5 


130 


1.63 


0±25.86 


3.86 • IO" 5 


2.36 


0±28.13 


2.26 • IO" 5 


140 


1.66 


0±25.72 


3.79 • 10" & 


2.46 


0±27.42 


2.12 • 10~ 5 


150 


1.61 


8.27±22.97 


4.25 • IO" 5 


2.49 


±24.85 


1.89 • 10~ 5 


160 


1.74 


0±14.91 


2.33 • IO" 5 


2.41 


0±20.96 


1.65 • IO" 5 


170 


1.48 


0±20.95 


3.46 • IO" 5 


2.35 


0±19.98 


1.61 • 10~ 5 


180 


1.42 


0±21.00 


3.60 • 10"° 


2.18 


15.0±13.35 


1.61 • 10" b 


190 


1.27 


0±16.19 


3.11 • IO" 5 


2.01 


0±29.29 


2.77 • IO" 5 


200 


0.85 


0±15.89 


4.54 • 10" b 


1.26 


6.49±15.96 


3.00 • 10" b 
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Table 5: Mass dependence of upper limits UL calculated using two runs together. 



M, MeV 


N bs ± cr 


UL 





0±11.21 


1.85 • 10~ b 


10 


0±11.05 


1.78 • 10~ b 


20 


0±10.56 


1.61 • 10~ b 


30 


0±9.72 


1.37 • 10~ b 


40 


0±8.62 


1.14 • 10~ b 


50 


0±7.42 


0.88 • 10~ b 


60 


0±6.39 


0.68 • 10~ b 


70 


0±5.78 


0.58 • 10- b 


80 


0±5.69 


0.49 • 10~ b 


90 


0±6.27 


0.50 • 10~ b 


100 


0±7.88 


0.58 • 10- b 


110 


0±12.59 


0.85 • 10~ b 


120 


0±25.67 


1.61 • 10~ b 


130 


0±19.91 


1.20 • 10~ b 


140 


0±15.51 


0.91 • 10~ b 


150 


0±14.94 


0.87 • 10~ b 


160 


0±16.79 


0.99 • 10~ b 


170 


0±17.07 


1.03 • 10~ b 


180 


0±20.09 


1.31 • 10~ b 


190 


0±22.12 


1.57 • 10~ b 


200 


13.66±19.06 


3.08 • 10~ b 
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Figure 1: K decay into sgoldstino and pions. 




Figure 2: Schematic view of the ISTRA+ setup. 
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Figure 3: Two photon effective mass spectra for events with one charged track and two showers in the SPi. 
Dashed line — first run data, solid line — second run data. 
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Figure 4: Missing energy spectra, second run data. Upper histogram — real data, middle histogram — back- 
ground MC events, lower histogram — signal MC events for mp = 90 MeV. 



10 



D 

Entries 

Mean 

RMS 



9908 
40356 
0.1370 



3000 



2500 - 



2000 



1500 - 



1000 - 



500 



1000 



800 



600 



400 



'00 



ID 9908 

Entries 12922 

Mean 0.1319 

RMS 0.1961E-01 



0.1 



0.12 0.14 0.16 



0.1 



0.12 0.14 0.16 



Figure 5: The influence of the "Veto" cut on the two-photon effective mass spectrum. Left part — second 
right part — first run; solid line — before the "Veto" cut, dashed line — after the "Veto" cut. 



11 




Figure 6: Missing mass squared distributions, real data, for successive cuts 3 — 10 respectively, second run data. 
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Figure 7: Missing mass squared spectra, background subtraction, first run — upper part, second run — lower 
part. Solid line shows spectrum for the neutral pion selection within narrow interval, dotted line shows spectrum 
for tails. Lower histogram in each part shows the difference. 
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Figure 8: Fit of the missing mass squared spectrum with sgoldstino signal at 120 MeV, first run data, upper 
part. Lower part — fit of the missing mass squared spectrum with sgoldstino signal at MeV, second run data. 
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Search for the decay K -> tt tt P 




Figure 9: Mass dependence of the upper limits, calculated using two runs together. 
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